Theoretical Search for Nested Quantum Hall Effect of Composite Fermions 
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Almost all quantum Hall effect to date can be under- 
stood as integral quantum Hall effect of appropriate particles, 
namely electrons or composite fermions. This paper investi- 
gates theoretically the feasibility of nested states of compos- 
ite fermions which would lead to a quantum Hall effect that 
cannot be understood as integral quantum Hall effect of com- 
posite fermions. The weak residual interaction between com- 
posite fermions will play a crucial role in the establishment 
of such quantum Hall states by opening a gap in a partially 
filled composite-fermion level. To treat the problem of in- 
teracting composite fermions, we develop a powerful method 
that allows us to obtain the low energy spectra at composite 
fermion fillings of v* — n + i> without making any assumption 
regarding the structure of composite fermions in the topmost 
partially filled level. The method is exact aside from neglect- 
ing the composite-fermion Landau level mixing, and enables 
us to study rather large systems, for example, 24 particles at 
a total flux of 62 hc/e, for which the dimension of the lowest 
Landau level Hilbert space is ~ 10 17 . We have investigated, 
for fully spin polarized composite fermions, several filling fac- 
tors between 1/3 and 2/5 using this approach. The results 
indicate that any possible incompressibility at these fractions 
is likely to have a fundamentally different origin than that 
considered earlier. 
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I. INTRODUCTION 

If electrons did not interact, only the integral quantum 
Hall effect (IQHE) [1] would occur in nature. The dis- 
covery of the fractional quantum Hall effect (FQHE) [2] 
signaled the existence of a new correlated state of mat- 
ter, the essence of which lies in the formation of quan- 
tum particles called composite fermions, [3] which are 
bound states of electrons and an even number of quan- 
tized vortices. Electrons avoid each other most efficiently 
by turning into composite fermions, which, in turn, in- 
teract much more weakly than electrons. For a large 
number of phenomena, it is a valid first approximation 
to neglect the interaction between composite fermions 
altogether, indicating that the interaction between com- 
posite fermions can often be treated perturbatively and 
does not cause any phase transitions. 

The fundamental property of composite fermions, 
which is responsible for the dramatic phenomenology of 
two-dimensional systems in high magnetic fields, is that 
as they move about, the Aharonov Bohm phase is partly 



canceled by the phase produced by the vortices tied to 
other composite fermions and as a result they experience 
an effective magnetic field B* given by 



B* = B- 2pp4>Q 



(1) 



where B is the external magnetic field, p is the electron 
density (in two dimensions), the even integer 2p is the 
vorticity of composite fermions (CF's), and 0o = hc/e. 
(The composite fermion with vorticity 2p is denoted by 
2p CF or CF-2p.) They form Landau-like levels in this re- 
duced magnetic field, with their filling factor (v*) related 
to the electron filling factor (y) by 



(2) 



2pv* ± 1 



The wave functions for composite fermions are given by 
[3] 



j<k 



(3) 



where z., = 



iyj denotes the position of the jth parti- 
cle, <5> u * is an antisymmetric wave function for fermions 
at the effective filling v* , and Plll is the lowest Landau 
level (LLL) projection operator. The wave function i Sf l> 
is of course the wave function of correlated electrons at 
v, but can also be interpreted as the wave function of 
composite fermions at filling factor v* , because on the 
right hand side, each electron has 2p vortices bound to it 
through the multiplicative factor TTj<j.(zj ~ z k) 2p , which 
converts each electron in into a composite fermion. 

Let us begin by neglecting the interaction between 
composite fermions. In this case, a gap opens up when 
the composite fermions fill an integral number of levels, 
i.e., when v* = n, which corresponds to electron filling 
factors given by 



(4) 



2pn+ 1 



where n is a positive or negative integer. [The negative 
values of n give v = \n\j {^2p\n\ — 1), and correspond to 
the situation in which the effective field B* is opposite to 
the external field B.\ The model of non-interacting com- 
posite fermions thus predicts quantum Hall effect (QHE) 
of electrons at these fractions, and also at 



1 



2pn + 1 



(5) 



1 



due to particle hole symmetry in the lowest Landau level. 
(The fractions in Eq. 5 can be obtained by formulating 
the problem in terms of holes rather than electrons in the 
lowest Landau level, and then making composite fermions 
out of them.) The QHE at these filling factors will be 
termed "CF-2p QHE." CF-2 QHE and CF-4 QHE are 
routine, and CF-6 QHE and CF-8 QHE have also been 
observed [4,5]. These sequences exhaust most of the more 
than 50 observed fractions to date in the lowest Landau 
level. The wave functions have been tested in detail for 
finite systems, and are very accurate, without involving 
any adjustable parameter [3,6-8]. 

In short, most of the QHE can be understood as the in- 
tegral QHE of appropriate fermions. It ought to be noted 
that the QHE of higher order composite fermions can be 
viewed as a consequence of interactions between lower 
order composite fermions. For example, the CF-4 QHE 
originates from interactions between CF-2's. Nonethe- 
less, the physics of all states at filling factors of Eqs. (4) 
and (5) is correctly described in terms of non-interacting 
composite fermions of the appropriate kind. 

While the above fractions are obtained most immedi- 
ately in the CF theory, it has been known since the very 
beginning of the composite fermion theory that other 
fractions are not ruled out. [3,9] New QHE between two 
successive integral Hall states of composite fermions will 
be called NQHE, where N stands for "nested." A simple 
scenario for NQHE is as follows. Consider electrons in 
the filling factor range 



1 



> v > 



2p{n + 1) + 1 ' 2pn + 1 
which map into 2p CF's in the range 
n + 1 > v* > n , 



(6) 



(7) 



containing n levels fully occupied and the (n + l) st level 
partially occupied. It is a natural conjecture that inter- 
CF interaction can possibly cause gaps at 



v = n + 



n 



2pn + 1 



(8) 



which will produce new fractions between the familiar 



fractions 



and 



rt+l 



Here, the 2p CF's in the 



2pn+l 2p(n+l) + l ' 

topmost partially filled level capture 2p vortices (as a re- 
sult of the residual 2p CF- 2p CF interaction) to turn into 
higher order composite fermions ( 2p+2p CF's) and con- 
dense into n levels, thereby opening a gap and produc- 
ing quantum Hall effect. This QHE state is a nested 
state with higher order composite fermions forming in 
the background of the original composite fermions. 

Certain qualitative consequences of the above scenario 
are as follows. The NQHE of composite fermions is anal- 
ogous to QHE of electrons between two successive IQHE 
states of electrons, that is, to FQHE in higher electronic 



Landau levels. It is known that FQHE is rare in second 
electronic Landau levels (i.e., in the range 4 > v > 2, 
including spin degeneracy) and non-existent in third or 
higher Landau levels; by analogy, NQHE is expected to 
be rare. Extending the analogy further, the strongest 
NQHE states for fully spin polarized composite fermions 
are expected to be at v* = 1 + 1/3 and v* = 1 + 2/3, 
which correspond to electron filling factors v = 4/11 and 
v = 5/13; these are expected to be the strongest nested 
CF states [3,9]. 

Recently, Pan et al. [10] have reported observation of 
4/11 and 5/13, heralding a new generation of quantum 
Hall effect that cannot be understood as the IQHE of 
composite fermions. This discovery has given fresh impe- 
tus to the issue of NQHE, especially because it is possible 
that more NQHE states will be observed in the future as 
the sample quality improves and the temperature drops, 
reminiscent of how the FQHE appeared on top of the 
IQHE. What other nested states are possible? What is 
the true nature of these states? Such questions have mo- 
tivated us to investigate this topic further. 

Even though the theoretical scenario described above 
is plausible, and indeed natural, there is no a priori guar- 
antee that it actually occurs, and it is important to carry 
out quantitative tests to ascertain its applicability to the 
real world. It is worth recalling here that not all frac- 
tions that can occur in principle do really occur in nature. 
For example, FQHE can occur for some hypothetical in- 
teraction at very small filling factors, but is preempted 
by Wigner crystal for the Coulomb interaction, and in 
higher Landau levels, the FQHE is often believed to be 
unstable to charge-density-wave type state. Therefore, it 
would require further theoretical work before one could 
claim with any degree of confidence that NQHE based 
on the physics described above is actually possible. 

The quantitative theoretical investigations so far sug- 
gest that there should be no FQHE at fractions other 
than those in Eqs. (4) and (5) for fully spin-polarized 
electrons. An exact diagonalization study [11] of N = 8 
electrons at 4/11 finds a non-uniform ground state, which 
is an evidence against FQHE liquid. In another approach 
[11,12], a model has been developed for the effective in- 
teraction between composite fermions in the second level. 
In Ref. [11] it was argued that no new FQHE is obtained 
[11]. Ref. [12] found that the form of this interaction 
favors a "bubble crystal" of composite fermions at fill- 
ing factors like 4/11 rather than the quantum Hall ef- 
fect. These studies indicate that the effective interaction 
between composite fermions is not sufficiently strongly 
repulsive to stabilize FQHE at 4/11 and several other 
fractions. Both approaches have their problems. In ex- 
act diagonalization, the system is effectively small, with 
only ~ 3 composite fermions in the second level, and one 
may question if the conclusion based on this system will 
hold up as N increases. The effective interaction model 
involves several assumptions, because of which the results 



2 



are plausible but not conclusive. 

In Ref. [13] showed that FQHE at 4/11 is very likely 
possible for a partially spin polarized system or a spin 
unpolarized system. It is not yet known for sure if the 
FQHE states at 4/11 and 5/13 are partially or fully po- 
larized in the experiments of Pan et at, [10] but given 
the relatively high magnetic fields at which these states 
have been observed, it is quite possible that they are fully 
spin polarized, and therefore it is worth revisiting the is- 
sue of whether fully polarized incompressible states can 
be found theoretically at these filling fractions. That is 
the prime motivation behind our present study. In this 
work, we will look for FQHE at several filling fractions 
in the range between 1/3 and 2/5. 

We will study finite systems of N particles at total 
flux 2Q (in units of the flux quantum 4>q = hc/e), where 
the relation between N and Q for a given filling factor 
will be derived from the physics described above. The 
most reliable theoretical method would be exact diago- 
nalization. However, because of the exponentially grow- 
ing Hilbert space, exact diagonalization is possible only 
for very small systems, and very few filling factors. For 
example, for v — 4/11, exact diagonalization has been 
possible for TV — 8 electrons. [14] The next system con- 
tains N — 12 (in the spherical geometry, discussed be- 
low), for which the lowest Landau level Hilbert space 
contains a total of 8.6 x 10 7 Slater determinant basis 
states; the diagonalization of matrices of such sizes is 
beyond the reach of present day computers. We have de- 
veloped a powerful new method, namely diagonalization 
in low-energy composite- fermion basis (LECFB), which 
we believe gives accurate results. It exploits the fact that 
the composite fermion theory allows us to directly iden- 
tify low-energy states of the full Hilbert space, which lets 
us work within a sub-space much smaller than the full 
Hilbert space. For example, we have been able to obtain 
the low-energy spectrum for the 4/11 state with as many 
as 24 electrons, where the dimension of the full Hilbert 
space is ~ 10 17 . 

A brief outline of the method is as follows. Consider 
the filling factor range given in Eq. (6). Given that the 
two ends of this range are well described as n + 1 and 
n filled levels of composite fermions, it is natural to ex- 
pect that the low energy states at v contain n composite 
fermion levels completely occupied, and the remaining 
composite fermions in the (n + l) st level. The configu- 
rations that involve promotion of composite fermions to 
higher kinetic energy levels are expected to cost substan- 
tially higher energy and are neglected. In other words, 
we include all configurations of type shown in Fig. (1) 
but neglect configurations of the type shown in Fig. (2). 
This neglect of composite-fermion Landau level mixing is 
the only assumption in our approach, which can be tested 
for consistency at the end of the calculation. Other than 
that the method is exact. The most important point 
here is that we do not make any assumption with re- 



gard to the interaction between composite fermions, or 
their structure in the partially filled level. In our calcu- 
lation, we first construct all states of composite fermions 
at v* — n + v, with the lowest n levels fully occupied 
and the next one with filling factor P, and then diag- 
onalizc the Coulomb Hamiltonian within this sector to 
obtain the spectrum of low energy states. The calcula- 
tion is rather involved, requiring extensive Monte Carlo, 
but still numerically stable, and yields reliable results. 

The present method is a significant advance over the 
variational method of Ref. [12]. It is much more reli- 
able because it eliminates several approximations made 
in that work. In Ref. [12], the problem of composite 
fermions at v* — n + v was mapped into that of fermions 
at v. Further, it was assumed that the interaction energy 
of the system of composite fermions in the partially filled 
level can be approximated by a sum of two-body interac- 
tions. The interaction between composite fermions was 
obtained by keeping only two composite fermions in the 
(n+l) st level and integrating out the composite fermions 
of the lower filled levels. In general, when there are many 
composite fermions in the (n + l) s * level, the integrating 
out of the lower-level composite fermions will produce a 
quite complex interaction, with two-, three-, and n-body 
terms, because the full system is strongly correlated. The 
hope was that the n-body terms do not cause any phase 
transitions, in which case the two-body terms will pro- 
duce the correct state. However, this approximation was 
untested. Finally, the method was a variational study, in 
which the energies of certain wave functions were com- 
pared to determine which had the lowest energy; there 
was of course no guarantee that this state would describe 
the true ground state. In contrast, in the present study, 
we work with the full composite fermion system and do 
not assume anything about the form of the interaction 
between composite fermions or the nature of their ground 
state. 

The following section contains a discussion of how we 
construct the low-energy composite fermion basis and 
how we carry out the Gram-Schmid orthogonalization 
and diagonalization. The reader not interested in the 
technical details can skip this section and directly pro- 
ceed to the subsequent section, wherein the results are 
given and their implications are discussed. 

II. DIAGONALIZATION IN LOW-ENERGY 
COMPOSITE-FERMION BASIS 

We first summarize certain relevant facts from the in- 
troduction. The electron filling factor v given by 



corresponds to composite fermion filling factor v* . We 
will be interested below in the specific values of v* given 
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by 



v = n + v = n + 



2pn+l 



(10) 



Here, composite fermions fill n levels completely and oc- 
cupy v fraction of the (n + l) st level. If the 2p CF's in 
the partially occupied level capture 2p additional vor- 
tices and condense into n filled Landau level state, then 
a gap would open up and NQHE would be obtained. The 
question is whether this mechanism really occurs in na- 
ture. 

In the composite fermion theory, there are several 
possible approaches for constructing a low energy ba- 
sis. In this work, we construct a LECFB of all states 
at v* = n + v in which n 2p CF levels are fully occupied 
and the next level has filling v. We do not make any 
approximation regarding the structure of 2p CF's in the 
partially filled level. This basis is in principle straight- 
forward: we simply need to take all electronic states at 
v* = n + v and attach 2p vortices to obtain the compos- 
ite fermion states. The method, however, is technically 
rather demanding, and the present section is devoted to 
it. 



A. Wave functions 

We will employ the spherical geometry in which N elec- 
trons move on the surface of a sphere under the influ- 
ence of a radial magnetic field B created by a magnetic 
monopole at the center. The magnitude of the B field is 
given by 2Qcf) /AiTR 2 where <f) — hc/e is known as the 
flux quantum, R is the radius of the sphere, and Q is 
called the monopole strength which should be either an 
integer or a half-integer because of Dirac's quantization 
condition. 

The interacting electron system at monopole strength 
Q maps into a system of weakly interacting composite 
fermions at an effective monopole strength q* = Q — 
p(N — 1). The wave functions for interacting electrons at 
Q are given by 



(11) 



where <E>i is the wave function for the fully occupied low- 
est Landau level, and (fr q * sue antisymmetric wave func- 
tions at q* . $ is in general a linear superposition of 
Slater-determinant basis states made up of the monopole 
harmonics, Y q * tS ^ m , given by [15]: 

Yq*,a,m{toj) = JV,.„ n (-l)«' + - m e i «'^ uf +m V f "™ 
2q* + s 



xV(-D" ; 



/ A 
r=0 



q* + s — m — r 



{v*V 3 y- r {u*U 3 Y , (12) 



N q* sm 



^(2q i_2^ + 1) (q 4 - 



m)\{q* 



m)\ 



1/2 



47T 



s!(2<7* + s)\ 



(13) 



s=0,l,2,... is the Landau level (LL) index, to be differ- 
entiated from n, the number of filled Landau levels. Qj 
represents the angular coordinates 6j and <frj of the jth 
electron, and 



Uj = cos(#j/2) exp(— i(j)j/2) 



vj = sin(6j/2)exp(i(j>j/2) . 



(14) 



(15) 



It was shown in Ref. [7] that 'J'q, the wave function of 
interacting electrons at Q — q* + p(N — 1), is obtained 
from <J> g * by replacing Yq*, s , m by , defined as: 

\rj \q* + s -m-rj 1 



_8_y (d_ 

duj) \dvj 



J! 



(16) 



where 



J 3 = Il^' Ufe ~ v i Uk ^ 



(17) 



Here the prime denotes the condition k ^ j. 



B. Low-energy composite- fermion basis 



For v = 0, that is at the special filling factors 



2pn± 1 



(18) 



the wave function for the ground state has the simple 
form: 



2pn±l 1 



2p,- 



(19) 



where $„ is the Slater determinant wave function for n 
filled Landau levels. Since $„ is unique, the wave func- 
tion ^ also contains no adjustable parameters. It 
has been tested in the past and was found to be remark- 
ably accurate [6,7]. 

Of interest in this work is the situation when v ^ 0. 
Here, we obtain the low energy spectrum through the 
following steps. 



where 
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1. We perform exact diagonalization for N particles 
at filling factor v in the lowest Landau level to get 
all eigenstate. The form of the interaction used for 
the diagonalization is of no significance, because its 
only role is to produce basis wave functions, but we 
work with the Coulomb interaction. Let us denote 
these states by 3>", where a labels different states. 

2. We then promote each eigenstate to the (n + l) s * 
Landau level, and fill the lower n Landau levels 
completely with additional particles to get the N 
particle state at v* = n + v. This gives all wave 
functions for which the lowest n Landau levels 
are fully occupied and the next Landau level has 
filling v. These states are denoted by L n $", where 
L n denotes addition of n Landau levels. It ought to 
be remembered that the operator L n changes the 
number of particles. 

3. We then multiply each wave function by <j> 2p and 
carry out the lowest Landau level projection by the 
method discussed above. This gives us correlated 
basis functions at v. These steps are summa- 
rized as 



n + v 



2pv* 



(20) 



(21) 



The advantage of constructing basis states in this 
way is that it directly gives us basis functions with 
well defined orbital angular momentum L, which is 
a good quantum number in the spherical geometry 
We can work within each L sector independently, 
because they are not coupled by the interaction. 
(Note that in the above equation, is a wave 
function for TV electrons only.) 

4. The basis thus obtained is not orthogonal. We 
obtain an orthogonal basis following the Gram- 
Schmid procedure. Some relevant details are out- 
lined below. 

5. We construct the Hamiltonian matrix. 

6. Diagonalization of the Hamiltonian gives the low 
energy eigenspectrum at v. 

If an incompressible state is obtained, the ground state 
wave function is likely to be well approximated by 



C. NQHE states on sphere 

Since our goal is to look the feasibility of NQHE de- 
scribed above, we do our calculations at those values of 
Q where it is possible for 2p CF's in the partially filled 
level to condense into a standard type of FQHE state. 
For that, we first need to determine the relation between 
Q and N for a given filling factor. Clearly, in the limit 
iV^oowe must have 



lim 

N — KX 



N 
2Q 



(23) 



where 2Q is the number of flux quantum. For finite sys- 
tems, however, the ratio N/2Q is not exactly equal to v 1 
and the identification of what filling factor a finite system 
corresponds to requires additional theoretical input. 

Let us first consider the incompressible states at filling 
factor v — 2 pii+i ' ^ nc wave function for which is given 
by Q^&n. For these states, the relation between Q and 
N can be obtained by noting that the product of two 
monopolc harmonics at q' and q" gives a monopole har- 
monic at q' + q", i.e. the monopole strengths add. For 
reasons that will become clear below, let us denote the 
number of particles by N and the monopole strength by 
q. The lowest LL projection is unimportant for the ques- 
tion of the relationship between N and q, as it does not 
alter q. In the spherical geometry, for monopole strength 
q the degeneracy of the lowest Landau level is 2q + 1 , for 
the next LL it is 2q + 3, and so on. From that, it is clear 
that the n filled LL state is obtained for q ft given by 



Qn 



N - 



2n 



In particular, one filled LL is obtained at 



N - 1 



9i = 



(24) 



(25) 



From the addition rule, the monopole strength of Q^&n 
is given by 

q = 2p Ql + q n =(p+l-^N-(p+^j (26) 

It can be verified that 2q/N gives the filling factor v in 
the thermodynamic limit. 

Now we ask what is the monopole strength q* for the 



state at v* = n + v = n + 



2pn+l 



for N electrons. Let 



us denote by N the number of fermions in the (n + 1) 
partially filled Landau level. In order for N fermions to 



PLLL$l P L n ^ s « PLLL$l P L n [P LLL $l P $ n ] (22) form the S3 StatC ' th ° de S enerac y in this level must 



where $p s « Plll^i^u is the wave function for the 
ground state for N particles at v = 2 pn+i ■ 



satisfy: 

2<f + 2n+l = 2g+l 
which implies that 



(27) 
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q =q-n 



(28) 



The total number of fermions, including the fermions in 
the lowest n fully occupied Landau levels is 



N = N + 2nq* + n 2 



(29) 



Eliminating N, one gets the following relation between 
N and q*: 



N = 2q* ( n + 



2pn + 1 / p + 



2 T +n 2 (30) 



In 



The state at v* = n + 2 pn+i * s obtained for all vaules of 
(N,q*) related by the above equation, provided TV is an 
integer and q* is an integer or a half integer. The state 
at v = 2 pv*+i m turn i s obtained for all values of (N,Q) 
with Q = q* +p(N -1). 

We have considered in this work the following filling 
factors. 

• v = Here, v* = 1 + g- The parameters are: 
p = p = n = n= 1, which gives 

» 3/V-8 n UN — 16 

<? = —5—. Q = — z — ( 31 ) 



• v 



j^: Here, v* = 1 + |. The parameters are: 



p = p = n = 1 and n = —2, which gives 

13JV- 17 

io 



* 3N ~ 7 n 

q = ^o~' 



• v = ^j: Here, v* = 1 + 



(32) 



Yg. Hue, ^ — j. -r g. The parameters are: 
p = p = n = 1 and n = 2, which gives 

* 5A " 17 , q= 19a :: 31 (33) 



14 



14 



• = j-: Here, z/* = 1 + g- The parameters are: 
p = n = n = 1 and p = 2, which gives 

,- = -^, Q = i-^- (34) 



12 



12 



Of course, only integer values are allowed for TV, 2q*, 
and 2Q. The table I explicitly lists the systems that we 
have studied below. 



D. Matrix elements 

The basis obtained in the way described above is not 
orthogonal for a given angular momentum sector L. We 
orthogonalize these using the Gram-Schmid procedure 
and than diagonalize the Coulomb Hamiltonian in the 
orthogonal basis. These require the knowledge of matrix 



elements of the type {4>i\4>ni) and (<fri\H\(j) m }, where <pi 
and 4> m (I to) are unorthonormalized wave functions, 
which we obtain by Monte Carlo. Since the Monte Carlo 
evaluation is most efficient when the integrand is positive 
definite, we determine these using the equations [16,17] 

^ = <M<fti + <hn) - {4>l\4>l) - {4>m\(j>m) 



and 



i\v\4> m ) = 



(35) 



+ <j> m \V\ct>l + <t>m) - {<t>l\V\<t>l) ~ {4m\V\<t>m) 



(36) 



where V denotes the Coulomb interaction Hamiltonian. 



E. Orthogonalization 

We denote {u{\ as the unorthogonal but normalized 
basis set obtained from the set {0;}, and define 



and 



Hj = ( u i\ u j) = 



(ui\V\uj) 



{4>i\4>o) 
\<l>i\\4>j\ 



(37) 



(38) 



\<kUi\ ' 
these are real numbers. 

We then follow the Gram-Schmid procedure to con- 
struct an orthogonal basis set {tpi} in terms of the un- 
orthogonal basis set {ui}. The former can be expressed 
as 

J- 1 j 

i=i * 
j'-i 

with Nf = (ipi\ipi) and fa = (ipi\uj)/Nf. While the 
above relation allows an iterative evaluation of {ipi}, we 
find it more convenient to work with the following ex- 
pressions: 

JV? = (^i) = l-E/M (40) 



fc=l 



and 



fij — 



N? 



1 

Nf 

1 

w 



i-1 



(«» - X! fk%1Pk\Uj) 
k=l 

i-1 



fe=l 



(41) 
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Eq. (39) was used to obtain these expressions, as well 
as the fact that Uij and are real. The Hamiltonian 
matrix elements 



Va = 
12 NiN-i 



(42) 



in the basis set are thus obtained as 
1 



3 NiNj 
i-i 



Vij - ^ fkjEikNiNk 
fe=i 

r i- 1 



- E /« { E U N i N i + E fkjEikNiNk 



i=i 



fc=i 



(43) 



which are again obtained iteratively. We finally diago- 
nalize the matrix to obtain the eigenvalue spectrum 
for a given L sector. 



F. Computation 

The major computation time is spent in calculating the 
matrix elements {<pi\cj) m ) and ((j)i\V\(j) m ) in the Monte 
Carlo method, because at each step we must evaluate 
M Slater determinants, M being the size of the LECFB. 
To give an example, at (N,Q) = (24,31) we have M — 
1656. To minimize the computational time, we sample all 
the states relative to the L = state, which is obtained 
from the lowest energy state at v* . We typically perform 
six to ten Monte Carlo runs for each (N, Q) with (0.6 - 
1.0) xlO 6 Monte Carlo iterations in each run. For larger 
systems, we use parallel computers. (We divide all of 
the Monte Carlo steps into several configurations with 
each configuration placed on a single node of a Bcowolf 
class PC cluster. One node consists of a dual 1 GHZ 
Pentium III Processor. To obtain one data point, at a 
particular N, we use as many as 30 nodes repeatedly until 
the standard deviation in energies is sufficiently low to 
produce the desired accuracy.) Approximately 80 CPU 
days were needed for the calculation of the spectrum for 
(N,Q) — (24,31). The computation time increases very 
quickly with N, approximately as iV 12 . 

It is difficult to ascertain how the statistical error in 
each element of the matrix would propagate into the 
eigenvalues. We instead obtain the spectrum in several 
different runs, and use the eigenvalues from those differ- 
ent runs to obtain the average and the standard deviation 
for each eigenvalue. 



III. RESULTS AND DISCUSSION 

We have obtained the low energy spectra at several 
values of (N, 2Q) which were identified earlier with sys- 
tems at special filling factors 4/11, 5/13, 7/19, and 6/17. 



Figs. (3), (4), (5), and (6) show results for 4/11; Figs. (7), 
and (8) for 5/13; Figs. (9), and (10) for 7/19; and 
Figs. (11), (12), and (13) for 6/17. The energy spec- 
trum for (8, 18) have been given in the past both from 
exact diagonalization and the composite fcrmion theory 
[6,8,11]. The accuracy of the energies is good enough to 
be able to state whether the ground state has L = or 
not. The error increases with increasing L, which is a 
consequence of our use of an L = state for sampling; 
in principle, the error at large L can be reduced by sam- 
pling with another reference state, but we do not see any 
need for it for the question of interest in this work. 

The orthogonalization makes a qualitative difference in 
the spectrum. Figs. (14) and (15) show the energies of the 
basis states prior to orthogonalization, which ought to be 
compared to the spectra in Figs. (5) and (6), respectively. 

The results explicitly confirm that our method in- 
deed gives very low energy states. For example, con- 
sider 4/11. The ground state energy is approximately 
— 0.419e 2 /do per particle, which compares very favor- 
ably with — 0.420e 2 /d 7 which is what one would get 
from a linear extrapolation between the energies [7] at 
1/3 (-0.410e 2 /d ) and 2/5 (-0.433e 2 /eZ )- Further, a 
large fraction of the states obtained here are within an 
energy band less than the effective cyclotron gap at 1/3, 
which is on the order of 0.1 e 2 /cIq. In other words, most 
of the states considered here lie within the 1/3 gap. This 
justifies our approximation of neglecting the states that 
involve transitions between composite fermion Landau 
levels, which are expected to cost energy of order 0.1 
e 2 /d . 

Now we ask whether the spectra imply incompressibil- 
ity. A characteristic feature of all incompressible states 
known to date is that they have an L = ground state, 
which is separated by a gap from other states; in addi- 
tion, there is also a well defined branch of neutral ex- 
citation. The ground state is extremely well described 
as a filled Landau level state of composite fermions, 
and the excitation as a particle hole pair excitation 
of composite fermion [6,8,18]. This is the case, with 
no exception, for for systems with (Q, N) related by 
Q = (1 ± l/2n)A - (1 ± ra/2), which are identified with 
v = n/(2n± 1). 

None of the spectra have these standard features of 
an incompressible state. Even in cases when the ground 
states has L = 0, there is no well defined CF-exciton 
branch. The final spectrum at v thus bears little resem- 
blance to the spectrum at v that we started out with, im- 
plying that the interaction between composite fermions 
in the partially filled level is quite different from that of 
electrons at v in the lowest Landau level. The scenario 
described in the introduction for NQHE is thus not borne 
out. 

The ground states at the flux values considered above 
are in general not incompressible uniform states with 
L = 0. Consider 4/11 for example. It has been known 
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[6,8,11] that for N — 8 particle system at Q = 9, which 
is identified above with 4/11, the ground state does not 
have L = 0. However, the system is effectively very small 
here, with only N — 3 composite fermions in the second 
level, and one may question if the conclusion will hold 
up as N increases. We find that for TV = 12 and 20 
the ground state has L = 0, but for N = 16 and 24 it 
does not. This provides an illustration for why a study 
of several values of N is crucial; if one had only studied 
N = 12, that would have led one to exactly the opposite 
conclusion. Another example is the 12 particle system at 
6/17, which has an L — ground state and might sug- 
gest FQHE at 6/17; however, the spectra for TV = 18 and 
24 at 6/17 explicitly clarify that this is only a finite size 
effect that does not survive in the thermodynamic limit. 

Our study implies that even though the ground state 
sometimes has L = 0, it does not describe an incompress- 
ible state at any of the filling factors considered here. The 
following observations support this statement: 

1. When the ground state has L ^ 0, the L = state 
is a rather highly excited state. 

2. In the plot of the ground state energy with 1/iV, 
the energies of the L = ground states are not 
particularly low. See Fig. (16) for 4/11. 

3. The ground state at v generally produces an ex- 
cited state at v, even in the L = sector (prior to 
orthogonalization), which indicates a breakdown of 
analogy between v and v. 

4. The energy per particle is not a smooth function of 
1/-/V, but instead shows fluctuations characteristic 
of a compressible state. See Fig. (16). 

5. It is not inconceivable for a compressible state to 
have L = 0, because a uniform state avoids Hartree 
energy cost, but it would be difficult to understand 
how an incompressible state might have L/0. 

As a result, the physical mechanism, in which some 
of the composite fermions turn into higher order com- 
posite fermions and condense into Landau levels to ex- 
hibit QHE, does not appear to be relevant for fully spin 
polarized electrons. The interaction between composite 
fermions in the second level is thus not of a form to pro- 
duce NQHE. It must be stressed, however, that the re- 
sults do not rule out NQHE through some hitherto un- 
known mechanism, which would imply a different relation 
between 2Q and N for a given filling factor than the one 
assumed above. 

In Ref. [12], which also investigated some of the above 
fractions, it was predicted that the ground state is a 
bubble crystal of composite fermions, similar to what 
was found in Hartree Fock studies of corresponding elec- 
tronic states in higher Landau levels [19]. Aside from 
several approximations, the results of that study were 



also predicated upon the quantitative validity of certain 
trial wave functions, making an independent confirma- 
tion worthwhile. Unfortunately, the spherical geometry 
used here is not convenient for the study of charge den- 
sity wave states. We do note in passing that a crucial 
step in Ref. [12] was to model the states at v = n + v 
in terms of fermions at v with an effective two-body in- 
teraction. A close similarity between our Fig. (11) and 
Fig. (5) of Ref. [11], obtained from the two-body interac- 
tion model [both for (A, 2Q) = (12, 30)] provides support 
to the validity of this model. 

In the end, we note that the above calculations assume 
a strictly two dimensional system (with zero width), no 
disorder, and no mixing with higher electron Landau lev- 
els. These approximations represent an idealization of 
the actual experimental system. It is expected that the 
idealized model, while not accurate quantitatively, should 
give the correct qualitative behavior. 
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TABLE I. The monopole strength Q as a function of TV 
for the NQHE states studied in this work. The quantities v 
and v* are the the filling factors of electrons and composite 
fermions, with v being the CF filling factor in the partially 
filled level. TV is the total number of composite fermions, and 
N is the number of composite fermions in the partially filled 
level, at filling V. Q and q* are the monopole strengths for 
the electrons and composite fermions. 
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FIG. 1. The configurations included in our study. Each 
composite fermions is depicted as an electron carrying two flux 
quanta. The lowest composite fermion level is fully occupied 
and the second one is partially occupied. 
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FIG. 3. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (12,29), identified with 4/11, where 
N is the number of composite fermions and 2Q is total flux 
penetrating the sample. The energies in this and subsequent 
figures are quoted in units of e 2 /elo, where e is the dielectric 
constant for the background material (e « 13 for GaAs) and 
lo = \J hc/eB is the magnetic length at the relevant filling 
factor (4/11 in this case). Spherical geometry is employed for 
the calculation, and L is the orbital angular momentum of 
the state. 
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FIG. 2. The configurations neglected in our study. These 
have higher "kinetic" energy than the configurations shown 
in Fig. (1). 
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FIG. 4. The low energy spectrum of interacting composite FIG. 6. The low energy spectrum of interacting composite 

fermions at (N, 2Q) = (16, 40), identified with 4/11. fermions at (N, 2Q) = (24, 62), identified with 4/11. 
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FIG. 5. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (20,51), identified with 4/11. 



FIG. 7. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (14,33), identified with 5/13. 
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FIG. 8. The low energy spectrum of interacting composite FIG. 10. The low energy spectrum of interacting composite 

fermions at (N, 2Q) = (19, 46), identified with 5/13. fermions at (N, 2Q) = (23, 58), identified with 7/19. 
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FIG. 9. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (16,39), identified with 7/19. 



FIG. 11. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (12,30), identified with 6/17. 
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FIG. 12. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (18,47), identified with 6/17. 



FIG. 14. The low energy spectrum of composite fermions 
at (N,2Q) = (20,51) prior to orthogonalization. The corre- 
sponding spectrum after orthogonalization is given in Fig. 5. 
The spectra shown here and in Fig. 15 are not physically 
meaningful, but are given here only to show the effect of or- 
thogonalization. 
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FIG. 13. The low energy spectrum of interacting composite 
fermions at (N,2Q) = (24,62), identified with 6/17. 
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FIG. 15. The low energy spectrum of composite fermions 
at (N,2Q) — (24,62) prior to orthogonalization. It should be 
compared to the spectrum in Fig. 6. 
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FIG. 16. The ground state energy per particle as a function 
of 1/N for 4/11. The energies plotted here have been mul- 
tiplied by the factor J Pn/p, where pjv is the density of the 
finite system and p is the density in the thermodynamic limit, 
to account for the dependence of density on N. It ought to be 
noted that the ground state does not have L — in general. 
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